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1  Introduction 


Optimal  assembly  of  a  coherent  system  from  a  given  set  of  n  components  is  an  important 
theme  of  research  in  reliability  Theory.  For  instance  see  El-Newcihi.  Proschan  and  Sethu- 
raman  (1986).  Derman,  Lieberman  and  Ross  (1974).  Boland.  El-Neweihi  and  Proschan 
(1988)  (henceforth  referred  to  as  BEP1).  Boland.  El-Neweihi  and  Proschan  (1992)  (hence¬ 
forth  referred  to  as  BEP2).  Shaked  and  Shantikumar  (1992)  (henceforth  referred  to  as  SS) 
among  others.  Such  an  assembly  may  involve  bolstering  the  original  components  of  the 
system  by  spares  which  can  be  used  as  replacements  upon  failure  ( standby  redundancy)  or 
connected  in  parallel  with  the  components  ( active  redundancy).  The  original  components 
can  also  be  strengthened  by  performing  a  certain  number  of  minimal  repairs  at  each  failure. 
BEP2  and  SS  have  demonstrated  the  relevance  of  various  modes  of  stochastic  ordering  of 
lifetimes  of  components  in  optimal  allocation  problems.  We  extend  their  results  to  more 
general  setups  using  mainly  the  familiar  techniques  of  AI  and  Schur  functions. 

BEP2  considered  the  optimal  allocation  of  a  single  spare  as  a  standby  redundancy  to 
one  of  the  n  components  of  a  series  or  parallel  system.  In  this  paper  we  study  the  optimal 
assignment  of  n  spares  which  are  available  as  standby  redundancy,  one  each,  to  the  n 
components  of  series  or  parallel  system,  thus  extending  the  results  of  BEP2.  These  results, 
given  in  Section  3,  may  be  described  informally  as  follows.  If  the  components  and  the 
spares  are  hazard  ( reverse  hazard)  rate  ordered  and  the  system  is  a  series  ( parallel)  system, 
then  one  should  assign  the  stronger  ( weaker )  spares  to  the  weaker  components,  in  order. 

We  also  consider  the  optimal  allocation  problem  in  other  systems.  In  one  class  of 
systems,  the  system-lifetimes  are  determined  by  damages,  introduced  by  shocks,  exceeding 
random  thresholds.  The  problem  here  is  the  optimal  matching  between  the  thresholds 
of  the  components  and  the  parameters  of  the  processes  governing  the  arrival  of  shocks. 
In  another  class  of  systems,  it  is  possible  to  bolster  the  system  by  performing  a  certain 
number  of  minimal  repairs  to  components.  When  the  available  number  of  minima]  repairs 
is  specified,  one  has  to  find  the  optimal  assignment  of  these  minimal  repairs  to  the  various 
components  to  maximize  reliability  of  the  system.  In  a  third  class  of  problems,  we  consider 
a  large  system  based  on  n  similar  modules.  It  is  possible  to  bolster  the  first  component  of 
each  module,  by  placing  spares,  in  active  redundancy.  Given  the  total  number  of  spares 
available,  we  find  the  optimal  allocation  of  spares  to  the  components.  These  optimal 
allocation  results  are  obtained  in  Section  4.  Some  of  the  results  here  are  new  and  others 
improve  upon  earlier  ones  obtained  by  SS  and  BEP1. 

We  collect  a  few  definitions  of  various  orderings  among  random  variables  and  their 
interrelationships  in  Section  2.  We  will  assume  that  the  reader  is  familiar  with  the  theory 
of  Schur  and  AI  functions  as  found  for  instance  in  the  textbook  of  Marshall  and  Olkin 
(1979). 
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2  Preliminaries 


We  will  now  define  some  of  the  standard  concepts  of  ordering  among  random  variables. 
In  this  paper  we  will  use  the  words  increasing  to  denote  nondecreasing  and  decreasing  to 
denote  nonincreasing.  Let  A’  and  Y  be  two  random  variables  with  distribution  functions 
(df's)  F  and  G. respectively.  Let  F{x)  —  1  —  F(x)  and  G(x)  —  1  —  G(.r).  We  say  that 

A  >  Y  if  F(t)  >  G(x)  for  all  x. 

hr  _ 

A"  >  V  if  F(x)/G(x)  is  increasing  for  all  x. 

rhr 

X  >  Y  if  F(x)/G(x)  is  increasing  for  all  x.  and 

It 

X  >  Y  if  f(x)/g[x)  is  increasing  for  all  .r 

where  /  and  g  are  the  probability  density  functions  (pdf's)  of  A  and  Y.  respectively.  It  is 
well  known,  for  instance  see  Ross  (1983).  that 

lr  hr  si 

X  >  Y  implies  that  A'  >  V  implies  that  A’  >  V,  and 

It  _  T'nr  _  _  st 

X  >  Y  implies  that  X  >  Y  implies  that  A  >  Y. 

When  X  and  Y  are  both  nonnegative  random  variables  with  absolutely  continuous 

hr-  _  rhr 

distributions,  then  saying  that  A'  >  Y  ( X  >  Y)  is  equivalent  to  saying  that  the  hazard 
(revers  hazard)  rate  function  of  X  is  pointwise  less  than  (greater  than)  or  equal  to  that 
of  Y .  This  explains  the  nomenclature  for  such  orderings.  Earlier  work  in  this  area,  for 
instance  see  Keilson  and  Sumita  (1982),  have  called  the  hazard  rate  ordering  as  uniform 
stochastic  order  in  the  positive  direction. 

Consider  a  unit  whose  lifetime  has  df  F.  The  concept  of  minimal  repair  upon  the  failure 
of  this  unit  that  has  been  successfully  used  in  Reliability  Theory  .  see  Barlow  and  Hunter 
(1960)  and  Ascher  and  Feingold  (1984).  and  can  be  described  as  follows.  If  the  item  fails 
at  time  t  then  the  minimal  repair  amounts  to  replacing  the  item  with  a  functioning  item 
of  the  same  age.  more  formally,  the  df  of  the  lifetime  of  the  repaired  item  is  given  by 
{F[. r)  —  F(t))/F(t)  for  x  >  t.  Let  M{t)  be  the  number  of  minimal  repairs  performed  in 
time  t.  if  minimal  repairs  are  performed  every  time  an  item  fails.  It  is  well  known  that 
{M(t),  t  >  0}  is  a  nonhomogeneous  Poisson  process  with  parameter  —  log  F(1 ),  see  Ascher 
and  Feingold  (1984). 
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3  Optimal  assignment  of  standby  redundancy  to 
components  of  series  and  parallel  systems 


We  begin  with  a  brief  summary  of  some  existing  work  on  optimal  allocation  of  redun¬ 
dancy  in  coherent  systems.  Consider  a  system  S  of  n  independent  components  with  lifetimes 

A’i,  A 2 ,  . . . ,  A'n.  Suppose  that  there  are  independent  spares  with  lifetimes  Vj.  V2 . Yn. 

The  components  can  be  enhanced  by  placing  the  spares  in  parallel,  one  to  each  component 
that  is  in  active  redundancy.  Suppose  that  the  components  and  spares  are  ordered  by  the 

st  st  st  st  st  st 

usual  stochastic  ordering,  that  is,  A’i  <  A'2  <  •  •  •  <  Xn  and  V’i  <  Y2  <  •  •  ■  <  Vr,  and  that  S 
is  a  A'-out-of-n  system.  BEP1  obtained  fairly  general  results  on  optimal  allocation  in  this 
problem.  They  showed  that  the  optimal  allocation  is  to  assign  spares  to  components  in  the 
reverse  order  to  maximize  the  reliability  of  the  system,  and  in  fact  that 

hk]n(P(XTr  V  V;;  >f),r  =  l,2,...,n) 

is  an  AD  function  of  (Tr.ir’)  .  where  h^n(pi,p2. . . .  ,pn)  is  the  reliability  function  of  a 
A--out-of-n  system  and  7r  and  7r"  are  two  permutations  of  (1,2, . . . ,  n). 

Instead  of  using  the  spares  in  active  redundancy,  one  can  use  them  as  replacements 
upon  failure,  that  is  in  T„andby  redundancy.  The  problem  of  optimal  allocation  of  spares 
in  this  mode  of  redundancy  has  proved  of  be  less  tractable.  In  the  recent  paper,  BEP2. 
the  following  setup  was  considered.  Suppose  that  there  is  a  single  spare  which  can  be 
assigned  to  some  component  as  a  replacement  upon  failure.  Suppose  that  the  components 

hr  hr  hr 

are  ordered  in  the  hazard  rate  ordering,  that  is  Ad  <  A'2  <  •  •  •  <  Xn  and  that  the  system  S 
is  a  series  system.  BEP2  showed  for  this  special  case  that  the  single  spare  should  be  assigned 
to  the  weakest  component.  For  the  case  of  a  parallel  system  with  components  ordered  by 
the  reverse  hazard  ordering,  they  also  showed  that  the  spare  should  be  assigned  to  the 
strongest  component  to  maximize  reliability.  They  also  provided  counterexamples  to  show 
that  the  hazard  rate  and  reverse  hazard  rate  orderings  cannot  be  weakened  to  stochastic 
ordering  and  that  their  results  do  not  extend  to  more  general  A--out-of-n  systems. 

In  this  section  we  consider  this  situation  where,  instead  of  just  one  spare,  there  are  n 
spares  to  be  assigned,  one  each,  in  standby  redundancy  to  the  n  components,  and  obtain 
optimal  allocation  results  for  both  parallel  and  series  systems. 

The  following  lemmas  concerning  random  variables  ordered  by  the  hazard  rate  and 
reverse  hazard  ordering  are  the  key  results  in  obtaining  our  optimal  allocation  results  of 
this  section. 


^  hr 

Lemma  1  Let  X ,Y  be  random  variables  with  df's  F.G  satisfying  X  >  Y' .  Let  b(x)  be  a 
bounded  nonnegaiive  increasing  function.  Then  for  any  bounded  increasing  function  h(: r) 


/  h.{x)b(x)dF(x)  f  h(x)b(x)dG(x) 
f  b(x)dF(x)  ~  f  b(x)dG[x) 
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Proof:  To  prove  the  (3.1)  it  is  enough  to  prove  that 

I  d~J  f  b(x)dG(x)  (  b{x)dF{x)  —  f  b(x)dF{x)  f  b(x)dG{x)  >  0 

for  all  t  and  for  the  class  of  functions  b(x)  of  the  form  b(x)  =  J2o<,<m  where  a,  >  0. 

-oc  =  to  <  ti  <  i2  <  ■  ■  ■  <  tm  <  tm+ 1  =  oc  and  It(x)  =  1  if  x  >  1  and  =  0  if  x  <  t.  Fix  a 
t  <E  (  — oc.oc)  such  that  tx  <  t  <  <t+1  for  some  i  with  0  <  ?  <  m.  In  the  following,  we  will 
observe  the  usual  conventions  that  summations  over  vacuous  regions  are  zero.  It  is  easy  to 
see  that 


J  =  E  aJ(G(iJ)-G(())][  ^  atF(4)+  E  QrF(tT) 

0<J<1  0< 

~  E  —  -^(0)  E  QkG(tk)  +  E  arG(tr) 


t+ 1  <r<m 


=  E  *?[ mGit^-GWFitj) 

0<j< i  1  J 

0 <}<k<i  L 

+  E  E  a,ar[F(<r)(G(iJ)-G(0)-G(fr)(JF(G)-F(0)l.  (3.2) 

0<j<ii+l<r<m 


Since  A'  >  Y,  we  get  the  two  inequalities 


F(t)G{tj)  ~  G{t)F{ti)  = 


m  G{t) 
F(t,)  G(t:) 


F{tj)G{tj)  >  0  when  j  <  ?  (3.3) 


Fiumt,)  -  G(t))  -  Giumt,)  -  F(t)) 

=  -F(tr)(G(t:)-G(t))  G(tr)(F(t,)-F(t)) 

Git,)  G(t j)  F{t j) 


F{tj)G{tj)  >  0  when  r  >  j. (3.4) 


Applying  the  inequality  (3.3)  to  the  first  two  terms  of  (3.2)  and  the  inequality  (3.4)  to  the 
last  term  of  (3.2),  we  find  that  I  >  0.  <0 

It  is  clear  that  (3.1)  holds  even  when  we  can  allow  the  functions  h  and  b  to  be  also 
unbounded  as  long  as  all  the  integrals  in  (3.1)  are  finite.  For  such  and  other  generalizations 
see  Caperaa  (1988). 

r  /}  f 

Lemma  2  Let  Ai,A2,Vi,V2  be  independent  random  variables  such  that  A'j  >  A’2  and 
Y\  >  Y2.  Then 


P( at  >  Y1)P( x2  >  y2)  >  P( A,  >  v2)P{ a2  >  v,). 
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Proof:  Let  the  df's  of  A'j,  A'2,  V,  V2  be  f\,  F2,  G\ .  G2.  respectively.  To  prove  (3.5)  we  see 
that 


p( x1  >  }\)P(x2  >  y2)  -  P(x,  >  \2)P(X2  >  \\) 

=  J  Gi(xi)dFi(x-i)  J  G2(x2)dF2(x2)  —  J  G2(xi)dF][x-i)  J  Gi(x2)dF2{x2)  >  0 

by  appealing  to  Lemma  1  and  taking  F  =  f\,G  =  F2.  b  =  G2  and  h  =  G\/G2.  0 

We  now  exploit  the  above  lemma  to  obtain  optimal  allocation  results  for  series  and 
parallel  systems  when  spares  have  to  be  allocated  to  components  as  standby  redundancy. 

Theorem  1  Consider  a  system  S  consisting  of  n  components  in  series  with  independent 

lifetimes  AT,  A'2, . . .  ,  Xn.  Suppose  that  n  independent  spares  with  lifetimes  )\A2 . V„. 

which  are  independent  of  the  components,  are  available  as  standby  redundancy,  one  each  to 

hr  _  hr  hr  hr  hr  hr  ^ 

a  component.  Suppose  that  A]  >  A2  >  •  •  •  >  An  and  Yy  >  V2  >  •  •  •  >  V„ .  For  t  >  0,  and 
for  permutations  ir.ir"  of  ( 1,2, ....  n),  define 

5i(7r'  7r*)  =  P({Xrr  +  V;.)  >  t.r=  1, . . . . n). 

Then  gt(- 7r,7r*)  is  an  AD  function  for  each  t  >  0. 

Proof:  By  standard  methods  for  AD  functions  it  is  enough  to  prove  this  result  for  the  case 

hr  rhr 

n  =  2.  Note  that  >'i  >  Y2  is  equivalent  to  —  V2  -j-  t  >  —Vi  +  t.  for  each  t.  From  (3.5)  it 
follows  that 


p( at  +  r2  >  t)P(x2  a  \\  >t)>  p(Xi  +  Fi  >  t)P(x2  +  y2  >  t). 

This  proves  that  <7((x,7r’)  is  an  AD  function  when  n  =  2.  0 

From  this  theorem  it  is  clear  that  in  order  to  maximize  the  reliability  of  the  system  S 
one  should  allocate  the  spares  with  lifetimes  Vi.  V2, . . . ,  Vn  to  the  components  with  lifetimes 
A'l .  A2. . . . ,  A’n  in  reverse  order,  that  is  the  stronger  spares  should  be  allocated  to  the  weaker 
components  in  order. 

Theorem  2  Consider  a  system  V  consisting  of  n  components  in  parallel  with  independent 
lifetimes  X\ ,  X2, . . . ,  Xn .  Suppose  that  n  independent  spares  with  lifetimes  Vj.  V2, . . . .  V„, 
which  are  independent  of  the  components,  are  available  as  standby  redundancy,  one  each  to 

rhr  rhr  rhr  rhr  rhr  rhr 

a  component.  Suppose  that  X\  >  X2  >  •  ■  ■  >  Xn  and  Y\  >  V2  >  •••  >  Yn.  Fort  >  0, 

and  for  permutations  7r,  7r"  of  (1, 2, . . . ,  n),  define 

gt(n,Tcm)  =  1  -  P({ ATr  +  Vr; )  <  f.r  =  1, - n). 

Then  gfn.ir’)  is  an  AI  function  for  each  t  >  0. 


5 


rhr 

Proof:  Again,  it  is  enough  to  prove  this  result  for  the  case  n  =  2.  Note  that  >  i’2  is 

hr 

equivalent  to  —  V'2  +  t  >  4- t.  From  (3.5)  it  follows  that 

P( Ai  4-  Yi  <  t)P(X 2  4-  Vi  <  <)  >  P{X\  4-  Vj  <  t)P(X2  4-  V2  <  *)■ 

This  proves  that  gt{Tr.Trm)  is  an  AI  function  when  n  =  2.  0 

From  this  theorem  it  is  clear  that  in  order  to  maximize  the  reliability  of  the  system  V 
one  should  allocate  the  spares  with  lifetimes  Vj,  V'2. . . . ,  )  „  to  the  components  with  lifetimes 

A'i,  A'2, _ X„  in  order,  that  is  the  stronger  spares  should  be  allocated  to  the  stronger 

components  in  order. 

Remark  1  In  Theorem  1  and  Theorem  2  if  we  take  Tj  to  be  any  nonnegative  random 
variable  and  take  Vjj  =  I3  =  •••  =  V„  =  0  we  obtain,  as  a  special  case,  the  results  of 
BEP2  for  the  optimal  allocation  of  a  single  spare  mentioned  earlier  at  the  beginning  of  this 
section.  BEP2  also  showed  that  if  the  hazard  (reverse  hazard)  rate  ordering  among  the 
lifetimes  of  the  components  in  the  series  system  S  ( V )  were  weakened  to  just  stochastic 
ordering,  then  such  an  optimal  allocation  result  will  not  be  true.  These  counterexamples 
also  show  that  we  cannot  weaken  the  conditions  of  hazard  (reverse  hazard)  rate  ordering 
in  Theorems  1  (Theorems  2). 

Remark  2  Suppose  that  the  series  (parallel)  structures  considered  in  Theorems  1  and 
Theorems  2  constituted  a  module  within  a  larger  coherent  structure.  The  optimal  allocation 
conclusions  of  these  theorems  continue  to  hold  for  the  larger  coherent  system  as  well. 

Theorems  1  and  Theorems  2  find  applications  when  we  know  that  the  lifetimes  of  com¬ 
ponents  and  spares  are  ordered  by  the  hazard  and  reverse  hazard  rate  orderings.  Such 
orderings  are  easily  apparent  if  the  lifetimes  belong  to  a  proportional  hazards  or  propor¬ 
tional  reverse  hazards  family.  For  another  family  of  distributions  see  the  following  corollary. 

Corollary  1  Let  U,V  be  two  random  variables  with  logconcave  density  functions.  Let 

r  =  (rj.  r2, . . . .  rn)  and  s  —  (s-[,S2, - sn).  Suppose  that  the  distribution  of  the  lifetimes  of 

the  independent  components  and  spares  are  given  by 

x.=  Y.  <4 •>:=  E 

1  =  1,2 where  the  random  variables  Uv  are  independent  copies  of  U  and  VtJ  are 
independent  copies  ofV.  Suppose  that,  the  spares  are  to  be  allocated  to  the  components,  one 
each,  as  standby  redundancy.  Then 

gt(r,s)  d=  P(( X,  -f  }')  ><),*  =  ! . n) 

is  an  AD  function  of  (r.s)  and 

ht(r .  s)  d=  1  —  P([X,  +  V')  <  t).i  =  1 - ,rc) 

is  an  A I  function  of  (r,s). 
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Proof:  It  is  well  known,  for  instance  see  Karlin  and  Proschan  (1960),  that,  if  U  has  a  log 

It  hr  "  rhr 

concave  density  then  A,  <  X3  whenever  r,  <  r:  and  hence  that  A’,  <  X3  and  A’,  <  Xj. 
A  similar  conclusion  holds  for  the  lifetimes  of  the  spares.  The  corollary  now  follows  from 
Theorems  1  and  Theorems  2. 

4  More  Applications  of  Schur  and  AI  Functions  in 
Optimal  Allocation 


In  this  section  we  study  some  further  optimal  allocation  problems.  Once  again  we  use 
the  standard  techniques  of  Schur  and  AI  functions  to  pinpoint  the  optimal  allocations.  It 
will  become  clear  that  these  techniques  constitute  the  most  natural,  and  at  the  same  time, 
very  powerful  tools  in  this  area. 

Consider  a  system  ST  consisting  of  n  components  situated  at  n  locations.  At  location 
i,  shocks  arrive  according  to  a  Poisson  process  {Ni(t),t  >  0}  with  parameter  A,  and  the 
j th  shock  produces  a  random  damage  Di:  to  the  component  situated  at  that  location,  and 
component  i  dies  when  the  total  accumulated  damage  becomes  greater  than  or  equal  to  its 
random  threshold  A’,,  ?  =  1, 2, ....  n.  In  other  words 

Ti  <t  if  £  Dij  >  X, 

where  T,  is  the  lifetime  of  component  i.  We  assume  that  the  random  variables  {DtJ  are 

i.i.d.  with  a  common  logconcave  density  and  that  {{A'i . A’„,  },  {Dij,  i  =  1.2,....  n,j  = 

1, 2  {N\(t).t  >  0}, . . . ,  {Nn(t),  t  >  0}}  are  independent.  We  will  call  such  a  system 

as  a  shock-threshold  system. 

The  following  theorem  pinpoints  the  optimal  matching  between  the  thresholds  and  the 
parameters  of  the  shock  processes  for  both  series  and  parallel  systems. 


Theorem  3  Consider  a  shock-threshold  system  ST  as  above. 

hr  hr 

(A )  Suppose  that  ST  is  a  series  system  based  on  its  components.  Suppose  that  A’i  <  AA  < 

hr 

•  •  •  <  An .  Let  Ts{\)  be  the  lifetime  of  the  system  ST.  Then  for  every  increasing  function 


E(g(Ts{ A))  is  AI  in  A. 


(B)  Suppose  that  ST  is  a  parallel  system  based  on  its  components. 


rhr 

Suppose  that  A'i  < 


rhr  rhr 

A'2  <  •  •  •  <  A’n-  Let  Tp(\)  be  the  lifetime  of  the  system  ST .  Then  for  every  increasing 
function  g. 


E{g(Tp( A))  is  AD  in  A. 


Proof:  Notice  that 


P(Ts(X)>t)  =  P(X,  >  J2  D'>'i  =  1* 

l<?<A\(f) 

=  Zp(x<>  L  A»,i  =  i 

where  k  =  (ki,k2, . . .  ,kn).  Now 

-P(A,  >  ^  =  1,  ...,n)  =  f...fp(  Xi  >  xt. - Xn>xn)1[[f{k‘){xi)dxi...dxn. 

To  say  that  the  A',’s  are  hazard  rate  ordered  is  the  same  as  saying  that  P( X,  >  x)  is  TP? 
in  i  and  x.  Since  f(x)  is  logconvave,  f^[x)  is  T  P2  in  k  and  A’,  from  Karlin  and  Proschan 
(1960).  From  Hollander,  Proschan  and  Sethuraman  (1977)  products  of  T  P2  functions  are 
AI  functions.  Hence  P( AN  >  z,,...,An  >  xn)  is  A1  in  x  and  I"[i >s  AI  in  ( x.k ). 
Thus  from  the  preservation  theorem  Hollander,  Proschan  and  Sethuraman  (1977)  (Theorem 

3.3)  P{  X,  >  DtJ.i  =  1 - ,  n)  is  AI  in  k.  Again  fl.  exp  —  A,t— A;-  is  AI  in  (X.k). 

By  invoking  the  same  preservation  theorem  we  obtain  that  P(Ts( A)  >  t)  is  AI  in  A.  This 
proves  part  (A)  pertaining  to  a  series  system.  Part  (B)  is  proved  in  a  similar  fashion.  {> 
Suppose  that  in  a  shock-threshold  system  described  above  we  could  permute  the  com¬ 
ponents  amongst  the  locations,  still  putting  only  one  component  at  each  location.  The 
above  theorem  shows  that  if  the  system  is  a  series  system  based  on  its  components  then  we 
should  allocate  the  weaker  components  to  the  locations  with  lesser  intensity  of  shocks,  in 
order,  to  stochastically  maximize  the  lifetime  of  the  system.  If  the  system  were  a  parallel 
system  based  on  its  components,  then  we  should  allocate  the  stronger  components  to  the 
locations  with  lesser  intensity  of  shocks,  in  order,  to  stochastically  maximize  the  lifetime 
of  the  system.  This  agrees  with  the  intuitive  notions  regarding  series  and  parallel  systems. 

We  will  now  consider  a  system  S  with  n  components  and  suppose  that  we  can  bol¬ 
ster  the  components  by  performing  minimal  repair  upon  failure.  SS  considered  this  prob¬ 
lem  when  the  lifetimes  of  the  components  are  ordered  stochastically  and  there  is  a  vector 
(k i,  k2. . . . .  kn)  representing  the  numbers  of  minimal  repairs  that  can  be  performed  on  the 
components.  When  <5  is  a  parallel  system  they  showed  that  the  optimal  allocation  is  to  as¬ 
sign  smaller  kfs  to  the  weaker  components,  in  order.  In  what  follows  we  provide  a  simpler 
proof  of  this  result  and  give  a  similar  resuK  when  S  is  a  series  system. 

St  S  St 

Theorem  4  Suppose  that  X\  <  X2  <,•••  <  Xn,  with  df's  Ft,F2. - Fn,  are  the  indepen¬ 

dent  lifetimes  of  components  1,2.  ..,n.  Let  k  =  (&],...,&„)  denote  the  number  of  minimal 
repairs  assigned  to  components  respectively.  Let  Hv(t.k)  (Hs(t.k)  be  the  survival 

distribution  of  a  parallel  (series)  system  formed  of  components  l,...,n.  Then 

(i)  ffp(t.k)  is  AI  in  k  for  each  t  >  0. 

(ii)  Hs(t.k)  is  AD  in  k  for  each  t  >  0. 


n )  II  exp  (  —  A,t) 


(A,/ 


k,\ 
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Proof:  From  t he  standard  techniques  used  to  prove  that  Junctions  are  AI  or  AD.  we  can 
assume  that  n  =  2.  Let  M,(t)  be  the  number  of  minimal  repairs  performed  in  time  t  on 
component  i.  i  =  1.2  assuming  that  we  can  perform  an  unlimited  number  of  repairs.  Then 

E(M\(t))  =  -log(l  -  ^i(O)  >  -log(l  -  F2(1))  =  E(U2(1). 

Let  A ,{k)  be  the  lifetime  of  the  ?th  component  when  k  minimal  repairs  are  performed  on 
it.  Let  A i  <  k2.  Since  the  Poisson  distribution  is  likelihood  ratio  ordered  in  its  parameter, 
it  is  also  hazard  and  reverse  hazard  rate  ordered  in  its  parameter.  Thus 

P(Miit)  >  Aq  +  1  )P{M2(t)  >  k2  +  1)<  P{Mx{t)  >  k2  +  l)P(M2(t)  >  A,  +  1). 


Now 


H-p{t,  Aq .  Aq>)  —  1  —  P( M\  ( t )  >  Aq  +  1 .  M2(k2)  >  k2  +  1 ) 

>  1  -  P(M}{t)  >  k2  +  l.M2(k2)  >  Aq  +  1) 

—  1  —  P(X\{k2)  <  t .  A;(/,'i )  <  t )  =  Hr(i.k2.  k} ). 

The  second  part  is  established  in  a  similar  fashion  by  using  the  fact  that  the  Poisson 
processes  M\(t)  and  M2(t)  are  reverse  hazard  ordered.  <N 

Remark  3  In  the  above  theorem,  if  we  were  given  a  vector  k  with  ki  <  k2  <  ■  ■  ■  <  kn 
and  we  could  choose  a  permutation  7r  and  allocate  kTi  minimal  repairs  to  components 

i,  i  =  1.2 . n.  The  optimal  allocation  that  maximizes  stochastically  the  lifetime  of  a 

parallel  (series)  system  based  on  these  components  is  to  choose  7r  =  (1.2 . n)  (7 r  = 

(n.n  -  1 . .  . . ,  1 )).  The  first  part  of  the  above  theorem  is  due  to  SS.  Our  proof  here  is 
shorter  and  our  use  of  A I  functions  allowed  us  to  extend  the  results  to  the  series  case. 

We  end  this  section  by  considering  a  situation  of  optimal  allocation  in  which  specific 
components  within  modules  of  a  system  are  targeted  for  bolstering  by  means  of  active 
redundancy.  We  describe  some  variants  of  this  problem  that  have  been  considered  before 
and  then  present  our  extension. 

Consider  rt  components  CY  C2 . Cn  which  can  be  bolstered  by  adding  spares  in  paral¬ 
lel.  that  is  in  active  redundancy.  Suppose  that  we  have  a  vector  (kj ,  k2 . kn)  representing 

the  number  of  spares  which  are  available  for  active  redundancy.  When  the  components  and 
spares  are  identically  distributed  and  the  system  is  a  A*-out-of-n  system,  BEPl  showed  that 

the  reliability  of  the  system  is  Schur-concave  in  (Aq.  k2 , - kn).  W'hen  these  components  are 

parts  of  7?  similar  -ubsvstems  and  the  larger  system  is  a  series  system  based  on  these  subsys¬ 
tems,  SS  showed  that  the  reliability  of  the  supersytsem  is  Schur-concave  in  (Aq.  k2 . kn). 

Such  results  immediately  pinpoint  the  optimal  choice  of  (Aq .  k2 . kn).  We  now  generalize 

both  these  results  to  the  case  where  the  components  are  parts  of  n  similar  subsystems  and 
the  larger  system  is  a  m-out-of-n  system  based  on  these  n  subsystems. 
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Theorem  5  Let  P  be  a  coherent  system  with  A’  components.  Let  the  first  component  of 

P  be  independent  of  the  rest  of  the  components.  Let  P\ .  P2 . Pn  be  ti  independent  copies 

of  P.  Let  C,  be  the  first  component  of  P,  and  suppose  that  we  wish  to  bolster  C,  with 
k\  independent  spares  in  active  redundancy,  i  =  1  , ...,n.  Suppose  that  the  spares  and 

the  components  C\ . Cn  have  lifetimes  with  common  nf  F.  Consider  a  larger  system  S 

which  is  an  m-out-of-n  system  based  on  P\.  P2 . Pr-  Let  T(k)  be  the  lifetime  of  S .  Then 

E(g(T(k)}  is  Schur  concave  in  k,  for  every  increasing  function  g. 

Proof:  Let  Tfik,)  be  the  lifetime  of  P,  when  its  first  component  C,  has  been  bolstered 
with  k,  spares  in  active  redundancy.  Note  that  P(T,(kt )  <  t)  =  a  -f  bFk,+l(t)  where  a.b 
are  nonnegative  constants  related  to  the  reliabilities  of  the  system  P  with  a  functioning  or 
nonfunctioning  component  1.  and  hence  do  not  depend  on  i.  Thus 

PiT[k)  ><)  =  !-  hm.n+l]n(a  +  bFk'+1(t).....a  +  bFkn+1{t)) 

where  6m_„+i|n(pi, . . .  ,pn)  is  the  reliability  function  of  an  m  —  7?  -f  l-out-of-r?  system  based 

independent  components  with  reliabilities  pj _ ,pn.  Pledger  and  Proschan  (1971)  showed 

that  . .  .  ,pn)  is  Schur  convex  in  (logpj - -  log  p„ )  and  is  clearly  an  increasing 

m 

function  of  its  arguments.  Let  k  >  k'.  Since  F(t)k+l  is  convex  in  k.  it  follows  that 
F(t)ki+l  +  F(t)k2+l  is  Schur  convex  in  (k\,k2).  We  can  use  this  fact  to  show  that  (log(o  + 

bFkl+1[t)),...,log(a  +  bFkn+1{t)))  >  (log(a  +  bFkl'+1{t)) . log  (a  +  bFkn+li{t))).  Hence 

P{T(k)  >  t)  <  P{T{k')  >  0  f°r  each  t.  Thus  E(g(T(k))  is  Schur  concave  in  k.  for  even- 
increasing  function  g 

The  above  theorem  generalizes  previous  results  of  BEP1  and  SS.  BEPl  obtained  the 
above  result  for  the  case  where  P  consisted  of  only  one  component.  For  more  general 
systems  P.  SS  obtained  the  above  result  when  the  larger  system  is  a  series  system  based 
on  Pi, ....  Pn,  that  is  when  m  =  n. 

Remark  4  Suppose  that  a  total  of  L  spares  is  available.  Since  E[g[T[k )  is  Schur  concave 
in  k  from  tbi  above  theorem,  we  have  to  distribute  the  spares  in  the  most  homogeneous 
fashion  to  obtain  the  optimal  allocation.  In  fact,  the  optimal  allocation  vector  is  given 
by  k  =  ( s  +  1 , . . . ,  s  +  1,5,...,  5)  where  s,r  are  determined  uniquely  by  the  equation  L  = 

r  n  —  T 

sn  +  r.  0  <  r  <  n  —  1 . 
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